Search for: All records

SUMMARY Elastodynamic Green’s functions are an essential ingredient in seismology as they form the connection between direct observations of seismic waves and the earthquake source. They are also fundamental to various seismological techniques including physics-based ground motion prediction and kinematic or dynamic source inversions. In regions with established 3-D models of the Earth’s elastic structure, such as southern California, 3-D Green’s functions can be computed using numerical simulations of seismic wave propagation. However, such simulations are computationally expensive, which poses challenges for real-time ground motion prediction and uncertainty quantification in source inversions. In this study, we address these challenges by using a reduced-order model (ROM) approach that enables the rapid evaluation of approximate Green’s functions. The ROM technique developed approximates three-component time-dependent surface velocity wavefields obtained from numerical simulations of seismic wave propagation. We apply our ROM approach to a 50 km $$\times$$ 40 km area in greater Los Angeles accounting for topography, site effects, 3-D subsurface velocity structure, and viscoelastic attenuation. The ROM constructed for this region enables rapid computation ($$\approx 0.0001$$ CPU hr) of complete, high-resolution (500 m spacing), 0.5 Hz surface velocity wavefields that are accurate for a shortest wavelength of 1.0 km for a single elementary moment tensor source. Using leave-one-out cross validation, we measure the accuracy of our Green’s functions for the CVM-S velocity model in both the time domain and frequency domain. Averaged across all sources, receivers, and time steps, the error in the rapid seismograms is less than 0.01 cm s−1. We demonstrate that the ROM can accurately and rapidly reproduce simulated seismograms for generalized moment tensor sources in our region, as well as kinematic sources by using a finite fault model of the 1987 $$M_\mathrm{ W}$$ 5.9 Whittier Narrows earthquake as an example. We envision that rapid, accurate Green’s functions from reduced-order modelling for complex 3-D seismic wave propagation simulations will be useful for constructing real-time ground motion synthetics and source inversions with high spatial resolution. 

Abstract The Turkana Depression, located between the Ethiopian and East African plateaus, displays an anomalous rift architecture. It is missing the narrow, magma‐rich morphology observed in the Main Ethiopian Rift that cuts through the Ethiopian Plateau. Instead, diffuse faulting and isolated volcanic centers are widespread over several hundred kilometers. Turkana has also experienced less magmatism over the last 30 Myr than adjacent plateaus, despite having a thin crust and residing above a mantle that is inferred to be hot and partially molten. We hypothesize that lithospheric weakening has been the key control on magma transport across the lithosphere in the Turkana Depression and subsequent rift development. Using poro‐viscoelastic–viscoplastic models of melt transport, we show that magma extraction across a thin, weakened lithosphere is slower than across a thick, elastic lithosphere. Our results suggest that pre‐rift lithospheric strength can explain the magma‐poor character of Turkana for most of its tectonic history. 

Abstract Deep earthquakes at depths below 500 km are under prohibitive pressure and temperature conditions for brittle failure. Individual events show diverse rupture behaviors and a coherent mechanism to explain their rupture nucleation, propagation, and characteristics has yet to be established. We systematically resolve the rupture processes of 40 large deep earthquakes from 1990 to 2023 and compare the rupture details to their local metastable olivine wedge (MOW) structures informed from thermo‐mechanical simulations in seven subduction zones. Our results suggest that these events likely initiate from metastable olivine transformations within the cold slab core and rupture beyond the MOW due to sustained weakening from molten rock at the rupture tip. Over half of the earthquakes likely rupture beyond the MOW boundary and are controlled by both mechanisms. Rupturing outside the MOW boundary leads to greater moment release, increased geometric complexity, and a reduction in rupture length, causing greater stress drops. 

Abstract Two‐phase flow, a system where Stokes flow and Darcy flow are coupled, is of great importance in the Earth's interior, such as in subduction zones, mid‐ocean ridges, and hotspots. However, it remains challenging to solve the two‐phase equations accurately in the zero‐porosity limit, for example, when melt is fully frozen below solidus temperature. Here we propose a new three‐field formulation of the two‐phase system, with solid velocity (v_s), total pressure (P_t), and fluid pressure (P_f) as unknowns, and present a robust finite‐element implementation, which can be used to solve problems in which domains of both zero porosity and non‐zero porosity are present. The reformulated equations include regularization to avoid singularities and exactly recover to the standard single‐phase incompressible Stokes problem at zero porosity. We verify the correctness of our implementation using the method of manufactured solutions and analytic solutions and demonstrate that we can obtain the expected convergence rates in both space and time. Example experiments, such as self‐compaction, falling block, and mid‐ocean ridge spreading show that this formulation can robustly resolve zero‐ and non‐zero‐porosity domains simultaneously, and can be used for a large range of applications in various geodynamic settings.